function d = hybridNeuron_AE(t, y, params)
    % hybridNeuron_AE: 神经元模型，第三方程含周期性电场 AE*cos(0.6*t)
    % y = [x; y; z]

    x     = y(1);
    y_var = y(2);
    z     = y(3);

    % 外部激励 us
    us = 0.37 * cos(0.7 * t);

    dxdt = (1 - params.a)*x - (1/3)*x^3 - y_var + us;
    dydt = params.c*(x - params.b*sin(params.beta*z)*y_var);
    dzdt = params.delta*y_var - params.alpha*z + params.AE*cos(0.6*t);

    d = [dxdt; dydt; dzdt];
end
